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RESUMEN 

El resumen sera traducido al espanol por los editores. We obtained nebular spectroscopy of LMC H n 
regions having classified stellar populations, thereby strongly constraining the ionization input parameters. 
Using photoionization models, we then evaluate the performance of nebular diagnostics of T* and abundance. 
We introduce [Ne m]/H/3 as a nebular diagnostic of the ionizing stellar T*. In contrast to the widely- used rf 
parameter, [Ne lli]/H/3 has greater sensitivity to mid and early O-stars, and is robust to nebular morphology 
and the presence of shocks. We present a preliminary calibration of both diagnostics for LMC metallicity. We 
also introduce 5*234 = ([S n] + [S in] + [S iv])/H/3 as a diagnostic of S abundance. 5*234 is much less sensitive 
to the nebular ionization parameter than is 523 or i?23. The intensity of [S iv]10.5/i is easily estimated from 
the optical and near-IR line ratios. We present calibrations of 523 and 5234 that are reliable at metallicitics 
Z < 0.5 Z Q . 

ABSTRACT 

We obtained nebular spectroscopy of LMC H n regions having classified stellar populations, thereby strongly 
constraining the ionization input parameters. Using photoionization models, we then evaluate the performance 
of nebular diagnostics of and abundance. We introduce [Ne iii]/H/3 as a nebular diagnostic of the ionizing 
stellar T+. In contrast to the widely-used rf parameter, [Ne ill]/H/3 has greater sensitivity to mid and early 
O-stars, and is robust to nebular morphology and the presence of shocks. We present a preliminary calibration 
of both T+ diagnostics for LMC metallicity. We also introduce 5234 = ([S n] + [S in] + [S Iv])/H/3 as a 
diagnostic of S abundance. 5234 is much less sensitive to the nebular ionization parameter than is 523 or R23. 
The intensity of [S iv]10.5/x is easily estimated from the optical and near-IR line ratios. We present calibrations 
of 523 and 5234 that are reliable at metallicities Z < 0.5 Z Q . 

Key Words: GALAXIES: ABUNDANCES — Hn REGIONS — STARS: FUNDAMENTAL PARAME- 
TERS — MAGELLANIC CLOUDS — SUPERNOVA REMNANTS 



1. INTRODUCTION 

Nebular emission-line diagnostics are one of our most important probes of conditions in cxtragalactic star- 
forming regions. In particular, these diagnostics are a primary means of constraining ionizing stellar populations 
and metallicities in distant galaxies. However, we rely heavily on photoionization models to interpret observed 
nebular emission-line ratios, and comparisons between theoretical models and observations have been lacking, 
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Fig. 1. Model photoionization tracks for rf vs [O m]/[0 II], using WR, 03, 07, and 09 stellar atmospheres. Symbols 
show our data, and numerals show nebular data from the Kennicutt et al. (2000) sample, with ionizing stellar spectral 
types binned according to the key; objects along the bottom are non-detections. See text for more details. 



especially for ordinary H n regions. Therefore, we have obtained long-slit spectroscopic observations of four 
H II regions in the Large Magellanic Cloud (LMC) , for the purpose of comparing the observed emission- line 
diagnostics with predictions from photoionization models. The sample spans a range in stellar spectral type 
from 07 to WNE, and in nebular morphology from classical Stromgren sphere to extreme shell. Two objects 
show evidence of shock excitation. 

We obtained both stationary, spatially-resolved observations of these nebulae, and scanned, spatially- 
integrated observations. A complete investigation of the entire dataset is presented by Oey et al. (2000; Paper I) 
and Oey & Shields (2000; Paper II). The spatially resolved observations generally showed excellent agreement 
between photoionization models and observations. However, a puzzling exception was our finding that the 
temperature-sensitive ratio [O iii]A4363/A5007 is owerpredicted by the models (Paper I). This is opposite to 
the expected effect of electron temperature fluctuations (Peimbert 1967). Our investigation of the abundance 
determinations shows that it is essential to use an appropriate relation between the electron temperatures in 
the high and low ionization zones, which could otherwise introduce errors in the abundance determinations of 
up to 0.2 dex. Interestingly, we find no evidence of abundance variations in DEM L199, a large Hn complex 
dominated by three early- type WR stars. For further details on our results from the spatially resolved data, 
we refer the reader to Papers I and II. Below, we summarize results from the spatially-integrated observations, 
which resemble those that would be obtained if the LMC were at a distance of 15-20 Mpc. 
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2.1. The rj' parameter 

The most widely-used diagnostic of T* has been the "radiation softness parameter" of Vflchez & Pagel 
(1988), defined as the relative ratios of singly to doubly-ionized O and S. 

,_ [OII]A3727/[OIII]AA4959,5007 
V = [S II] A6724/ [S III] AA9069, 9532 ' ( ' 

Figure [j] shows theoretical tracks of 77' vs. [O m]/[0 n], with our spatially integrated observations overplotted. 
The models are generated with the photoionization code Mappings II (Sutherland & Dopita 1993), using 
stellar atmospheres from Schaerer & de Koter (1997) and Schmutz etal. (1992). The default models (solid 
lines) have LMC metallicity and central hole size of 0.1 times the Stromgren radius (Rs)- The three tracks 
have electron density varying between 1, 10, and 100 cm~ 3 , which are directly proportional to variations in the 
ionization parameter U, with increasing U corresponding to increasing log [O iii]/[0 n]. The individual models 
represent spectral types of roughly 09, 07, 03, and early WR, toward decreasing values of 77', respectively. 
We also show model tracks for a shell morphology having an inner radius of 0.9 R s (dotted lines), and tracks 
having SMC metallicity but otherwise the same as the default. The data points for our sample of LMC nebulae 
are plotted with symbols identified in the key. Numerals show objects from the Kennicutt et al. (2000) sample, 
with large and small numerals indicating LMC and Galactic objects, respectively; these data are not spatially 
integrated over the entire nebulae (see Paper I). 

For both the models and data, Figure [l] shows that 77' can effectively discern the presence of WR stars 
and late-type O stars. However, mid- and early-type O stars yield highly degenerate values of 77'. Thus, this 
diagnostic is essentially insensitive to > 40, 000 K, with the exception of the WR stars. In addition, Figure [j] 
shows that 77' is also sensitive to nebular morphology. The model tracks for a hollow shell strucure (dotted 
lines) are clearly offset from the default models. 

Pcimbert etal. (1991) showed that the inclusion of supernova remnants (SNRs) within integrated H 11 
region spectra can affect interpretation of the nebular conditions. The object DEM 243 encompasses s SNR, 
and Figure |l| shows data points for the H 11 region both with and without the SNR included (open and solid 
diamonds, respectively). DEM 301 (cross) also shows evidence of shock excitation (Paper I). It is apparent in 
Figure |] that the 77' values of DEM 243 with the SNR, and DEM 301, overestimate T* in view of the actual 
stellar spectral types. This therefore shows that 77' is also sensitive to the presence of shocks. 

2.2. The [Ne m]/H@ parameter 

To address these shortcomings in the 77' diagnostic, we therefore introduce [Ne m]/H/3 as a complementary 
diagnostic of T+. Although this ratio is sensitive to abundance, the high ionization potential (40.96 eV) for 
Ne in yields greater sensitivity to hotter ionizing spectral types. Likewise, this high ionization potential reduces 
sensitivity to shock excitation. 

In Figure j2 we show model tracks of log [Ne lll]/H/3 vs log [O m]/[0 11], for sequences in T*, with symbols 
as in Figure [l]? Figure || shows that [Ne m]/H/3 is a fairly successful discriminant of T*, even between mid and 
early O-types. The dotted lines also show that it is highly insensitive to nebular morphology. In addition, the 
two data points for DEM 243 including and excluding the SNR (solid and open diamonds, respectively) now 
show agreement in [Ne m]/H/3, demonstrating insensitivity to shocks. This is further supported by the data 
for DEM 301 (cross), which also falls in the locus expected for its T*. Thus, [Ne m]/H/3 holds promise as a 
useful T* diagnostic for hot O stars. Interestingly, the Galactic and LMC objects do not seem differentiated 
as expected for their different metallicities. Table [l] gives a preliminary empirical calibration for these T+ 
diagnostics. Further details can be found in Paper I. 

3. DIAGNOSTICS OF METALLICITY: 5234 

H11 region metallicities are often estimated using the "bright-line" diagnostic of O/H: R23 = ([O 11] + 
[O m])/H/3 (Pagel etal. 1979). Recently, a similar diagnostic has been introduced for S: 5*23 = ([S 11] + 
[S m])/H/3 (Vflchez & Esteban 1996). Empirical calibrations for 523 have been presented by Christensen 
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Fig. 2. Same as Figure [j] for [Ne m]/H/3 vs [O Hl]/[0 n] 
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etal. (1997) and Diaz & Perez-Montero (2000), who suggested that it would be less [/-sensitive than i?23 and 
would also present a larger dynamic range. 

Figure || shows theoretical tracks for the abundance diagnostics, using the same photoionization code and 
stellar atmosphere models as before. These are overplotted with nebular data points from Dennefeld & Stasihska 
(1983), who obtained optical and near-IR spectrophotometry of a large sample of nebulae, with electron tem- 
peratures derived from the ratio A4363/A5007. Our models show that £23 is actually more sensitive to U 
than is R23. This results from the omission of [S iv], which has roughly the same ionization potential (35 
eV) as [O in]. The plot in Figure ||a, of log(S/H) vs. log 523, shows the large spread in the three tracks. 
These correspond to log U = —2, —3, and -4, with the first having the lowest values of S*23, as expected for a 
significant population of [S iv]. Note that this sequence of tracks is opposite to that for R23 (Figure ||c), which 
more completely samples the important ions. 

We therefore suggest an abundance diagnostic that includes [S iv], to better sample the relevant ions, and 



TABLE 1 

EMPIRICAL CALIBRATION OF T* DIAGNOSTICS 



Sp. type 


log [Ne m]A3869/H/3 


log 7?' 


WR 


> -0.6 


< 0.2 


03 - 04 


-0.9 to -0.6 




05 - 08 


-1.5 to -0.9 




09 and later 


< -1.5 


> 1.0 



a For LMC metallicity. 



NEBULAR EMISSION-LINE DIAGNOSTICS 




Fig. 3. Theoretical tracks for a) S23, b) 5234, c) 7223, and d) log (O/H) vs. log 5234, computed from models at 0.05, 
0.1, 0.3, 0.5, 1.0, and 2.0 Zq. Dashed, solid, and dotted lines correspond to U — —2, —3, and -4, respectively, and 
are overplotted with data from Dennefeld & Stasihska (1983). Straight lines are calibrations: solid line is fitted to our 
models for Z < 0.5 Z ; dashed is from Diaz & Perez-Montero (2000); dash-dot from Christensen etal. (1997). 
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thereby improve the effectiveness of the S diagnostic (Paper II): 

S234 = ([SII]A6724 + [S III] AA9069, 9532 + [S IV]10.5^)/H/3 (2) 

Although [S iv]10.5/x is a mid-IR line, its intensity can be estimated from a simple relation between optical 
and near-IR lines, for metallicities Z < 0.5 Zq (Paper II): 

[SIVll0.5u , [OIII]AA4959,5007 

l0g [S III] AA9069 9532 = ~°- 984 + [OII]A3727 < 3 > 

This relation provides a simple way to estimate S234. Figure [|& shows S234 for the same dataset, with the 
intensity of [S iv]10.5yLt computed from equation ^. The theoretical tracks are dramatically less sensitive to U 
than for S23, or even i?23 (Figure ||c). 

For S'234, a fit to our models gives a theoretical calibration (Figure |^>): 

log (S/H) = -5.58 + 1.27 log S'234 , (4) 

For S23, we have a similar theoretical calibration from the models (Figure ||a, straight solid line): 

log(S/H) = -5.43+ 1.33 logS23 . (5) 

The dashed and dot-dashed lines in Figure ||a show calibrations from Christensen etal. (1997) and Diaz & 
Perez-Montero (2000) for comparison. Figure ||d shows that caution is necessary in inferring O abundances 
from the S diagnostics, since apparently there is significant scatter in S/O. Figure || also shows that all of these 
calibrations are reliable only for Z < 0.5 Z . 

We thank Angeles Diaz, Rob Kennicutt, and Fabio Brcsolin for providing access to their data in advance of 
publication. MSO gratefully acknowledges financial support from the conference organizers that enabled her 
attendance at this celebration. 
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